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1.  INTRODUCTION  AND  GENERAL  NOTIONS 


The  purpose  of  this  report  is  to  provide  a  brief  review  of  the 
principles  governing  the  phenomenon  commonly  called  "velocity-coupled 
combustion  instability."  Most  of  the  basic  ideas  covered  here  origi¬ 
nated  many  years  ago.1  3  They  have  been  clarified  by  subsequent  work 
and  the  analytical  framework  into  which  they  fit  has  developed  consi¬ 
derably  in  the  recent  past. 

So  far  as  practical  applications  are  concerned,  the  general  problem 
of  velocity  coupling  remains  largely  unsolved.  The  most  difficult  ob¬ 
stacles  are  associated  with  the  experimental  work  required  to  produce 
the  values  of  certain  characteristic  quantities  arising  in  the  analysis. 
No  consideration  will  be  given  to  those  matters  or  to  the  details  of  the 
combustion  processes  and  origins  of  the  physical  behavior. 

Velocity-coupled  combustion  instability  is  most  simply  described  as 
the  behavior  under  unsteady  conditions  corresponding  to  erosive  burning 
in  steady  flow.  Observations  of  burning  solids  have  established  that  in 
steady  state,  the  regression  rate  is  dependent  on  both  the  pressure  of 
the  environment  and  on  the  rate  of  flow  parallel  to  the  surface.  For 
the  purposes  here  it  is  adequate  to  assume  that  the  burning  increases 
if  either  the  pressure  or  the  parallel  flow,  or  both,  are  increased,  as 
sketched  in  Figure  1.  The  symbol  Vii  denotes  the  velocity  parallel  to 
the  surface;  often  the  mass  flow  is  used  as  the  independent  variable, 
but  the  velocity  will  serve  well  enough  here  to  develop  the  principles. 


FIGURE  1.  Influences  of  Pressure  and  Parallel  Velocity 
on  the  Linear  Burning  Rate. 
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Now  imagine  that,  from  some  initial  steady  operating  point,  the 
pressure  is  slowly  changed  by  an  amount  6p.  If  r  —  p  ,  then  to  first 
order  in  small  quantities,  the  corresponding  change  in  burn  rate  is 
given  by 


6r 

r 


n 


6p 

P 


(1) 


where  r,  p,  are  the  initial  steady  values.  Obviously  6p  can  be  positive 
or  negative,  and  if  the  pressure  slowly  oscillates  (i.e.,  with  very  low 
frequency  or  long  period)  about  its  average  value,  then  so  does  the 
burning  rate.  Similar  behavior  occurs  if  the  parallel  velocity  is 
changed.  Figure  2  shows  schematically  these  quasi-steady  responses  of 
the  regression  rate  to  oscillations  of  pressure  and  velocity.  Figure  3 
shows  the  pressure,  velocity,  and  burning  rate  as  functions  of  time  for 
the  two  cases.  Note  that  for  slow  quasi-steady  oscillations,  the  burn¬ 
ing  rate  is  always  in  phase  with  whatever  independent  variable  is 
changed,  providing  the  burning  rate  increases  monotonically  with  the 
variable  under  steady  conditions.  We  restrict  ourselves  here  to  the 
simplest  situation  in  which  only  one  variable  changes  and  when  the  velo¬ 
city  fluctuation  is  smaller  than  the  average  velocity. 


FIGUKK  ...  Quasi -Steady  Changes  of  the  Une.tr  Burning  Hale. 


h 


FIGURE  3.  Quasi-Steady  Oscillations 
of  the  Linear  Burning  Rate. 


Consider  now  the  situation  in  a  combustion  chamber.  Suppose  that 
for  some  reason  the  chamber  pressure  p  is  changed  and  that  the  pro¬ 
pellant  surfaces  are  so  oriented  (e.g.  as  in  an  end  burner)  that  fluc¬ 
tuations  of  the  parallel  velocity  are  negligible.  When  the  pressure 
increases,  so  does  the  burning  rate,  which  causes  the  total  mass  of  com¬ 
bustion  products  in  the  chamber  to  be  increased.  If  no  other  process 
occurred,  the  pressure  would  be  further  increased,  an  unstable  situa¬ 
tion.  Of  course  in  practice,  the  increased  chamber  pressure  causes  the 
exhaust  flow  to  increase,  compensating  for  the  increased  burning  rate; 
it  is  then  a  familiar  result  that  if  n  <  1,  the  flow  in  the  chamber  is 
stable  to  fluctuations  of  pressure  when  the  exhaust  nozzle  is  choked. 
This  is  the  quasi-steady  limit  of  "I*"  or  "bulk  mode"  instability. 

Similar  considerations  apply  to  the  case  when  the  chamber  pressure 
is  held  constant,  but  the  parallel  velocity  changes.  The  condition  for 
stability  is  then  not  so  simply  stated  and  will  not  be  worked  out  here. 
However,  it  is  clear  that  for  the  cases  sketched  above,  an  increase  of 
parallel  velocity  will  cause  an  increase  in  the  rate  of  mass  addition  to 
the  chamber.  In  the  absence  of  other  processes,  this  will  cause  an  in¬ 
crease  of  chamber  pressure.  The  assumed  constraint  of  constant  chamber 
pressure  will  therefore  be  satisfied  only  if  other  processes,  including 
the  exhaust  flow,  are  accounted  for. 

In  any  practical  situation,  changes  of  the  velocity  parallel  to  the 
surface  will  not  occur  alone:  they  will  he  accompanied  by  changes  of 
pressure  as  well.  Moreover,  the  geometry  of  combustion  chambers  is  such 
that  a  quasi-steady  increase  of  pressure,  uniform  in  the  chamber,  will 
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produce  a  quasi-steady  increase  of  parallel  velocity  everywhere.  The 
latter  will  of  course  not  generally  be  uniform  over  the  burning  surface. 
However,  for  small  fluctuations,  the  total  change  of  burning  rate  will 
be  the  sum  of  the  changes  due  to  the  fluctuations  of  pressure  and  velo¬ 
city.  This  behavior  is  sketched  for  one  position  on  the  burning  surface 
in  Figure  4. 


FIGURE  4.  Quasi-Steady  Oscillations  in 
Response*  to  Both  Pressure  and  Velocity. 
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Qualitative  examination  of  quasi-steady  behavior  is  evidently 
straightforward,  although  the  details  of  a  particular  case  may  become 
tedious.  Under  truly  unsteady  conditions,  at  higher  frequencies,  the 
analysis  becomes  considerably  more  involved.  However,  one  principle*  is 
always  true:  an  increase  of  local  burning  rate  with  an  increase  of  local 
pressure  is  a  destabilizing  influence.  Thus  it  is  that  part  of  the 
burning  rate  in  phase  with  the  fluctuations  of  pressure  which  may  pro¬ 
duce  unstable  motions  in  the  chamber.  This  is  true  whatever  may  be  the 
cause  for  the  burning  rate  to  change. 

All  arguments  concerning  the  influence  of  velocity  coupling  on  com¬ 
bustion  instability  have  been  based  on  the  principle  stated  in  the  pre¬ 
ceding  paragraph.  What  sometimes  obscures  the  issue  is  the  necessity  to 
introduce  some  representation  of  the  way  in  which  the  unsteady  burning 
rate  is  connected  to  velocity  fluctuations.  This  problem  will  be 
treated  later  in  greater  detail  but  it  may  be  helpful  here  to  outline 
some  of  the  main  joints.  We  shall  consider  the  simplest  and  most  impor¬ 
tant  case  of  longitudinal  (axial)  acoustic  modes  in  a  straight  chamber. 
Even  with  combustion,  particulate  matter  in  the  gases,  average  flow, 
exhaust  nozzle,  and  other  features  of  a  real  motor,  the  unsteady  motions 
usually  found  may  be  regarded  in  good  approximation  as  acoustical  mo¬ 
tions  in  a  chamber  having  the  same  shape,  enclosed  by  a  rigid  boundary 
excluding  the  nozzle  volume  and  without  all  perturbations  associated 
with  combustion.  Thus  we  may  consider  the  lateral  burning  surface  to  be 
subject  to  classical  acoustic  modes  appropriate  to  a  closed/closed  tube; 
the  pressure  field  for  a  standing  wave  is  proportional  to  p  cos(knz) 

cos(w  t)  where  k  =  n/l/L  and  u)  =  ak  (n  =  0, 1 ,2, .  .  . ) .  To  be  definite, 

n  n  n  n 

we  restrict  our  attention  to  the  fundamental  or  first  axial  mode  having 
amplitude  c,  the  maximum  value  of  p’/p  : 


p'(z,t)  =  r.p  cos(Hj-)  cos  (jj~  t) 


(2) 


The  momentum  equation  for  acoustical  motions  is 


(3) 


It  will  be  shown  in  Section  2.2  that  this  principle  is  closely 
related  to  "Rayleigh’s  Criterion",  .»  statement  concerning  the  influence 
of  heat  .ild  it  ion  on  the  stability  o!  acoustic  waves. 
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and  with  equation  (1) 

-  3u’  -  7i  .  ,7i  .  ,nl  . 

p  at"  =  £p  l  •“  (E  z)  cos  (r t} 


This  is  easily  integrated  to  give 


u’(z,t)  =  e  i  sin  (j  z)  sin  t)  (4) 


which  can  be  written  to  show  explicitly  the  phase  difference  between  the 
pressure  and  velocity: 

u'(z,t)  =  e  i  sin  (j  z)  cos  (5) 

Comparison  of  equations  (2)  and  (5)  shows  that  for  z  <  L/2,  the  velocity 
lags  the  pressure  in  time  by  n/2.  To  appreciate  these  formulas  it  is 
helpful  first  to  sketch  the  pressure  and  velocity  mode  shapes  at  quarter 
cycle  intervals  (Figure  5).  It  is  also  useful  to  have  the  pressure  and 
velocity  fluctuations  as  functions  of  time  in  the  two  regions  z  <  L/2 
and  z  >  L/2  which  may  be  regarded  as  the  fore  and  aft  sections  of  a 
rocket  motor.  These  are  sketched  in  Figure  6. 

It  is  apparent  that  the  velocity  lags  the  pressure  by  n/2  in  the 
head-end  (0  <  z  <  L/2)  of  the  chamber  and  leads  the  pressure  by  n/2  in 
the  aft-end.  The  reason  why  the  phase  difference  depends  on  position  in 
the  chamber  may  be  seen  from  Figure  5.  For  the  fundamental  mode,  the 
velocity  is  always  zero  at  the  ends  only,  while  the  pressure  field  has  a 
node  at  the  center.  Thus,  although  the  velocity  at  any  given  instant  is 
in  the  same  direction  throughout  the  chamber-- i .e. ,  has  the  same  sign-- 
the  pressure  changes  sign.  That  sign  change  is  reflected  in  the  phase 
difference  sketched  above. 

Now  let  us  consider  the  response  of  a  burning  surface  to  such  an 
acoustic  field.  We  shall  base  the  discursion  on  several  assumptions: 

(1)  The  burning  rate  is  proportional  to  the  pressure  fluctuation 
for  pressure  coupling; 

(2)  The  burning  rate  is  proportional  to  the  velocity  fluctuation 
for  velocity  coupling; 
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FIGURE  6.  Pressure  and  Velocity  Oscillations 
in  the  Fore  and  Aft  Regions  of  a  Tube. 


(3)  Pressure  and  velocity  coupling  are  independent  and  linearly 
additive; 

(4)  The  proportionality  constants  for  pressure  and  velocity  coup¬ 
ling  are,  in  general,  functions  of  frequency  and  contain  the 
fact  that  the  change  of  burning  rate  is  not  in  phase  with  the 
pressure  and  velocity  fluctuations.  For  the  present,  rectifi¬ 
cation  is  ignored,  but  will  be  discussed  later. 
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There  ere  several  equivalent  ways  to  write  the  oscillations  of 
burning  rate  (items  (1)  and  (2)).  The  most  efficient  representation 
involves  complex  notation.  For  the  pressure  field,  write 


El 

P 


p  iu»t 
e 


(6) 


^  =  £  COS  (j— ) 


(7) 


It  is  necessary  to  choose  a  definition  of  zero  phase  and  maintain  that 
definition  consistently  throughout.  Unless  otherwise  specified,  we 
shall  follow  the  generally  accepted  convention  that  all  phases  are  mea¬ 
sured  with  respect  to  the  pressure  oscillation  so  p/p  is  always  real. 
The  actual  (observable)  values  are  found  by  taking  real  parts,  for  then 
equation  (2)  follows  from  equations  (6)  and  (7).  The  velocity  is: 


u  ’  u  iu)t 
—  =  -  e 


(8) 


.  a)  i?  .  (i) 

r-‘<  f  ■  *(•)  ■  • 2  (!) 


(9)  a 


where 


-  ( i  ) 


(!)  -1 


(9 )  b 


Note  that  because  the  acoustic  velocity  lags  the  pressure  by  nf 2,  u/a 
is  imaginary.  Define  the  response  functions  for  pressure  and  velocity 
coupling : 


Rb  =  Rb(r)  +  iRbU)  =  lRblei,t,b  =  lRbl  cos\  +  *b  (l0) 


V  V 


(i)  = 

(i)  . 

i*l>, 

l*vl' 
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Then  the  changes  of  burning  rate  due  to  pressure  and  velocity  coupling 
are  represented  as 


|Rb|  £  ei(u.t+^) 

P 


(12) 


.  .  (i)  i(u)t+i|J  + 

W  "  '  "  ’ 


(13) 


The  real  parts  of  these  are: 


/r'\(r)  6 

/-)  =  |Rb|§  cos  (u.t  +  ^ 

\  r  /p  p 


|R^|  cos  £  cos  u)t  -  |R^|  sin  i|^  ^  sin  wt 
P  P 


(14) 


r’\(r)  /uv(i)  n 

^(v)  =  cos(lut  +  2  +  V 


-  (i) 

=  |Rvl/“]  [cos(u)t  +  t[)  cos  4>v  -  sin(u)t  +  ^)sin  t^] 


=  lR 


q  (i) 

I  “cos  4>v  sin  wt  -  sin  cos  wt] 
a  1 


(15) 


The  par*  *  of  equations  (14)  and  (15)  which  are  proportional  to 
cos  wt  and  theretore  can  be  in  phase  with  the  pressure  oscillation  are: 


(16) 


(17) 


12 


NWC  TP  6363 


Note  that  according  to  equation  (9),  (u/a)^  has  a  negative  sign, 
(u/a)^  =  -(e/y)  sin(nz/L). 

What  these  results  show  is  that  for  pressure  coupling,  it  is  the  real 
part  of  the  response  function  which  controls  the  destabilizing  tendency 
and  for  velocity  coupling  it  is  the  imaginary  part. 

It  must  be  emphasized  that  the  phase  of  the  pressure-coupled  re¬ 
sponse  function  is  that  for  the  oscillatory  burning  rate  with  respect  to 
the  pressure,  but  the  phase  t|j  of  the  velocity-coupled  response  function 
is  for  the  burning  rate  with  respect  to  the  velocity  fluctuation.  The 
reason  that  the  imaginary  part  of  R  is  required  to  specify  the  de¬ 
stabilizing  tendency  of  velocity  coupling  is  that  the  acoustics  intro¬ 
duces  a  n/2  phase  difference  between  the  velocity  and  pressure  oscilla¬ 
tions,  as  shown  by  equation  (5).  Thus,  an  additional  phase  difference, 
provided  by  the  response  function  R  ,  is  required  to  produce  a  component 
of  the  burning  rate  which  may  be  in  phase  with  the  pressure.  Only  in 
this  way  is  the  principle  cited  earlier  satisfied. 

Choosing  the  parts  of  (r'/r)  and  (r’/r)  which  vary  as  cos  u)t 
serves  to  eliminate  those  portions ^which  cannotVaffect  the  stability  of 
pressure  waves.  According  to  the  principle  stated  on  page  7,  the  fluc¬ 
tuations  of  burning  rate  must  be  in  phase  with  p/p  to  be  destabiliz¬ 
ing  tendencies.  This  is  a  statement  about  the  sign  of  r'  ,  namely  that 
^  C  r ) 

r  /p  must  be  positive  for  instability.  Thus,  in  order  that  the  con¬ 
tributions  represented  by  equations  (16)  and  (17)  cause  instability: 


(r/r) ^ 

=  Rh(r)  > 

P/P 


(18) 


w;>v'’  ■(.)  .  , 

v 


(19) 


p/p 


P/P 


Note  that  for  quasi-steady  beha . ior,  the  sign  of  (r/r)  /(p/p)  is  fixed 
by  the  slope  of  the  steady  burn  rate  graphed  as  a  function  of  pre:,si:re 
or  parallel  velocity.  For  the  examples  shown  in  Figures  2  and  3,  the 
contributions  are  destabilizing. 

For  pressure  coupling  to  be  an  unstable  influence,  the  general 
statement,  according  to  equation  (18),  is  that  the  real  part  of  the 
pressure-coupled  response  function  must  be  positive.  This  is  true  for 
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all  frequencies  for  real  propellants.  If  R^  >  0,  all  elements  of 
burning  surface  everywhere  in  the  chamber  contribute  destabilizing  in¬ 
fluences.  For  velocity  coupling,  with  the  four  assumptions  listed 
above,  the  product  of  the  imaginary  part  of  the  velocity-coupled  re¬ 
sponse  times  the  ratio  of  the  spacial  parts  of  the  velocity  and  pressure 
fluctuations  must  be  positive  to  cause  instability.  The  situation  is 
more  complicated  than  for  pressure  coupling,  and  two  cases  must  be  con¬ 
sidered,  distinguished  by  the  sign  of  R 


R 


v 


(i) 


>  0 


For  velocity  coupling  to  be  destabilizing,  the  velocity  fluctuation 
must  have  a  component  in  phase  with  the  pressure  fluctuation,  u/p  >  0. 
Consider  again  the  fundamental  longitudinal  mode,  and  Figure  5;  u/p  is 
positive  only  in  the  head  end  of  the  chamber,  0  <  z  <  L/2.  Thus,  if  the 
imaginary  part  of  the  velocity  coupled  response  function  is  positive, 
the  contribution  from  burning  surface  in  the  forward  half  is  destabiliz¬ 
ing  and  the  contribution  from  the  aft  half  is  stabilizing.  This  conclu¬ 
sion  follows  directly  from  equation  (19)  upon  substituting  equations  (7) 
and  (8),  giving 


-  R 


(i)  (u/a_)(i) 
V  (P/P) 


»R<‘> 


TtZ 

sin  r 

nz 

y  cos  j~~ 


i  »y>  u„  (f) 


The  tangent  function  is  positive  in  the  forward  half  of  the  chamber, 
(i) 


R 


0 


Obviously  the  situation  is  the  reverse  of  that  for  R 


(i) 


>  0.  The 


contribution  from  the  forward  half  s  L  a  b  i  ^_i  n  g  and  the  contribution 
from  tlie  aft  half  is  destabilizing,  for  the  fundamental  mode. 

The  same  sort  of  reasoning  can  he  applied  to  the  higher  modes  of 
t fit*  chamber.  Suppose  that  the  imaginary  part  of  the  velocity  coupled 

response  is  positive,  ^  f  0.  The  pressure  fluctuation  always  has 

maximum  amplitude,  and  the  velocity  fluctuation  is  zero,  exactly  at  the 
head-end,  for  all  modes.  lit  a  region  near  the  head-end,  the  velocity 
lags  the  pressure  by  n/2.  Thus,  there  is  a  length  of  propellant  sur¬ 
face,  starting  at  the  head-end,  which  under  the  influence  of  velocity 
coupling,  is  destabi 1 izing  for  all  inodes. 
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At  the  aft  end  of  the  chamber,  the  pressure  fluctuation  also  has 
maximum  amplitude;  it  is  out  of  phase  with  the  pressure  at  the  head  end 
for  even  modes  and  in  phase  for  odd  modes.  The  velocity  fluctuation  at 
the  aft  *nd  leads  the  pressure  by  n/2  for  some  distance  upstream  of 
the  end.  There  is,  consequently,  some  portion  of  propellant  surface, 
stretching  from  the  aft  end,  which  is  stabilizing  under  the  influence  of 
velocity  coupling.  Regions  of  stabilizing  and  destabilizing  influence 
alternate  along  the  chamber,  and  can  readily  be  identified  by  extension 
of  the  preceding  argument.  Figure  7  shows  the  situation  for  the  third 

and  fourth  modes,  with  R  >  0.  To  compress  the  information  into 
single  figures,  the  mode  shapes  of  pressure  and  velocity,  are  shown  as 
they  occur  for  their  respective  maximum  amplitudes,  i.e.,  they  are 
sketched  for  instants  of  time  separated  by  one-quarter  of  a  period  (see 
Figure  5) . 


_  p(z) 

_ u(z) 

/  /  /  REGIONS  FOR  DESTABILIZATION  BY  VELOCITY  COUPLING  IF  Rv(i*  >  0 

FIGURE  7.  The  Influence  of  Velocity  Coupling  for 
the  Third  and  Fourth  Longitudinal  Modes. 


So  tar  in  this  discussion  nothing  has  been  assumed  concerning  the 

In  practice,  it 


relationship  between  the  response  functions  R.  and  R 

b  v 


is  often  supposed  that  R 


This  approximation  is 


v  is  proportional  to  R^. 

based  on  the  ideas  that  processes  in  the  gases  respond  in  a  quas i -steady 
fashion,  if  the  frequencies  of  oscillation  are  not  Loo  high,  and  that 
the  principal  source  of  a  time  lag  in  the  response  is  unsteady  heat 
transfer  in  the  solid  phase.  11  we  further  assume  that  the  response 
**b  ^°r  Rresiiure  coupling  is  reasonably  approximated  by  the 
two  parameter  form,4  then  we  are  able  to  relate  the  sign  of 

Figure  8  is  a  sketch  of 


function 
f ami  liar 

R  (i) 

V 

the 


to 
rea  1 


the  frequency  of  the  mode  considered, 
and  imaginary  parts  of  R.  .  and  therefore  R 


imaginary  parts  of  R^f  and 
the  dimensionless  frequency  0  =  wK/r^. 


as  functions  of 
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FIGURE  8.  Real  and  Imaginary  Parts 
of  the  Response  Function. 


The  real  part  R^ 
(i) 


(r) 


has  a  maximum  at  Q  =  at  which  frequency, 


or  nearly  there,  R,  passes  through  zero.  For  frequencies  less  than 
-2  D  (i) 

r  ^/k,  the  imaginary  part  of  Rb,  and  therefore  Ry  as  well,  is  posi¬ 
tive;  this  is  the  first  case  discussed  above.  Combined  with  the  discus¬ 
sion  of  the  first  and  second  cases  above,  this  completes  the  basis  for 
the  following  statements  often  encountered,  in  one  form  or  other,  in  the 
literature  dealing  with  velocity  coupling: 

-2 

(1)  In  the  low  frequency  range  ( i . e .  ,  u>  <  r  U^/k,  the  burning  rate 

leads  the  pressure  (i.e.,  Rb^*\  >  0)  and  velocity  coup¬ 

ling  in  the  forward  half  of  the  chamber  is  destabilizing  for 
the  fundamental  mode. 


-2 

(2)  In  the  high  frequency  range  (w  >  r  Qp/x)  the  burning  rate  lags 
the  pressure  (R^*\  R^*^  >  0)  and  velocity  coupling  in  the  aft 
half  of  the  chamber  is  destabilizing  for  the  fundamental  mode. 


An  important  aspect  of  velocity  coupling  which  has  so  jar  been 
ignored  in  this  discussion  is  rectification.  The  burning  propellant 
should  normally  be  sensitive  to  the  magnitude  of  the  parallel  velo¬ 
city,  but  not  the  direction.  Suppose  that  an  element  of  surface  is 
exposed  only  to  an  oscillatory  component  of  velocity  having  frequency  u>. 
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u’  =  u  cos  u)t.  Then  because  the  combustion  processes  should  be  sensi¬ 
tive  to  the  magnitude  but  not  the  direction  of  parallel  flow,  the  burn- 
ing  rate  should  be  related  to  the  magnitude  of  u': 


|  u1 1  =  | U  |  COS  UJt 


(20) 


which  is  sketched  in  Figure  9. 


FIGURE  9.  Rectification  of  a  Velocity  Oscillation. 


The  oscillatory  part  of  |u'|,  and  therefore  the  burning  rate  re¬ 
sponding  to  the  parallel  velocity,  has  frequency  2u),  not  u>.  Conse¬ 
quently,  the  input  disturbance  of  frequency  u>  cannot  be  reinforced  and 
velocity  coupling  cannot  be  acting  in  this  case.  This  of  course  is  due 
to  the  rectification  noted  above  and  shown  in  Figure  9.  Note  that 
rectification  also  produces  a  DC  component  which  may  be  reflected  as  a 
change  in  steady  erosive  burning. 

The  phenomenon  of  rectification  must  always  be  present,  but  it  does 
not  preclude  the  possibility  of  having  velocity  coupling  at  the  funda¬ 
mental  frequency  u>.  What  is  required  is  that  the  mean  velocity  in  the 
same  direction  as  the  fluctuation^  be  non-zero.  The  extreme  rase, 
sketched  in  Figure  10,  occurs  when  u  >  u'  ;  there  is  then  obviously  a 

component  of  |u  ♦  u*  |  at  the  frequency  w,  namely  u*  itself. 
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FIGURE  10.  Rectification  Prevented  by  Mean  Flow. 

If  0  <  u  <  |u'|,  the  situation  is  somewhat  more  complicated,  as 
shown  by  Figure  11.  The  part  of  |u  +  u'|  which  has  the  freqeuncy  u>  of 
u'  must  be  found  by  Fourier  decomposition  of  |u  +  ul.  Because  the  re¬ 
sults  depend  on  the  amplitude  of  u' ,  rectification  introduces  nonlinear 
behavior.  Note  that  this  is  a  consequence  of  the  kinematics  and  the 
assumption  that  the  response  of  the  burning  surface  is  sensitive  to  the 
direction  but  not  the  magnitude  of  the  velocity;  this  has  nothing  what¬ 
ever  to  do  with  the  details  of  the  combustion  processes. 

Now  the  average  velocity  u  varies  along  the  chamber,  in  general 
increasing  with  distance  from  the  head  end,  so  that  the  amplitude  of 
that  part  of  |u  ♦  u'|  having  frequency  u>  also  varies  along  the  chamber. 
The  earlier  discussion  remains  valid  if  the  definition  of  R  is  assumed 
to  include  the  influence  of  rectification.  That  is,  there  is  an  implied 
mult ipi icat ive  factor  proportional  to  the  Fourier  coefficient  associated 

with  that  component  of  I u  ♦  u’l  having  the  frequency  of  the  mode  being 

considered.2  3  r*  This  has  the  important  consequence  that  the  function 
Ry  is  a  strong  function  of  position  along  the  axis  of  the  chamber,  gen¬ 
erally  having  greatest  magnitude  in  the  aft  portion  of  the  chamber  where 
the  mean  flow  speed  is  larger.  Thus  one  should  anticipate  that  the  in¬ 
fluence  of  velocity  coupling  is  in  some  sense  more  sensitive  to  changes 
in  the  aft  end  of  the  chamber. 
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FIGURE  11.  Partial  Rectification  in  the  Presence 
of  Mean  Flow. 


For  example,  assume  that  the  dependence  of  the  velocity-coupled 
response  function  on  frequency  is  that  shown  in  Figure  8.  Suppose  that 
the  burning  rate  and  geometry  are  such  that  the  dimensionless  frequency 
ft  is  less  than  ft  .  The  influence  of  velocity  coupling  in  the  forward 

half  of  the  chamber  is  destabilizing  for  the  first  mode.  But  this  is 
weak  because  of  the  weakness  of  rectification  over  the  forward  portion 
of  the  grain.  Now  suppose  that  the  burning  rate  is  decreased  so  that 
ft  >  ft  ,  so  that  velocity  coupling  in  the  aft  end  where  rectification  is 

stronger,  is  deslabi 1 izing.  Such  a  change  of  burning  rate  may  be  accom¬ 
panied  by  a  sufficient  increase  in  the  influence  of  velocity  coupling, 
as  to  cause  au  initially  stable  fundamental  mode  to  become  unstable. 

These  remarks  serve  to  emphasize  that  the  influence  of  velocity 
coupling  is  inlr ins ii al ly  dependent  on  position  in  the  chamber.  It  is 
true  that  the  most  basic  piece  of  information  is  the  imaginary  part  of 
the  velocity-coupled  response  function,  such  as  sketched  in  Figure  8. 
This  should  ideally  be  known  as  a  function  of  dimens  ion less  frequency, 
without  thr  influence  of  rectification.  To  obtain  a  quantitative  as¬ 
sessment  of  the  stability  of  wave  motions  in  a  chamber  with  velocity 
coupling,  a  much  more  elaborate  analytical  structure  is  required  to 
account  far  spacial  variations  in  the  chamber.  That  is  the  suhjecl  of 
the  following  section. 
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2.  STABILITY  ANALYSIS  WITH  VELOCITY  COUPLING 


In  the  preceding  section,  the  answer  to  the  question  of  stability 
rested  ultimately  on  the  principal  stated  earlier.  However  reasonable 
that  proposition  may  seem — and  it  should  be  at  least  appealing,  because 
it  is  true--nevertheless  it  was  simply  stated  without  proof.  The  formal 
analysis  outlined  in  the  following  paragraphs  will  produce  an  explicit 
statement  of  the  principle,  and  in  such  a  form  that  it  can  be  used  to 
obtain  quantitative  results  for  practical  situations.  It  is  not  the 
purpose  here  to  provide  a  detailed  exposition  of  the  analysis,  but 
rather  to  review  the  salient  features.  First  the  linear  stability 
analysis  will  be  covered  and  then,  for  future  reference,  an  approximate 
nonlinear  analysis  will  be  summarized. 


2.1  FORMULATION  OF  LINEAR  STABILITY  ANALYSIS 

The  most  recent  treatment  of  linear  stability  analysis,  essentially 
the  form  used  in  all  current  computer  codes,  appeared  as  Reference  6.  A 
treatment  of  velocity  coupling,  using  essentially  the  same  analytical 
framework,  appeared  earlier  in  Reference  7.  The  idea  on  which  the 
analyses  are  based  is  that,  for  linear  behavior,  any  arbitrary  distur¬ 
bance  in  a  chamber  may  be  represented  by  Fourier  analysis  as  a  super¬ 
position  of  contributions  from  all  of  the  normal  modes  of  the  chamber. 
If  one  of  the  modes  should  be  unstable,  i.e.,  grow  in  time,  then  so  is 
the  disturbance.  Hence,  the  stability  of  an  arbitrary  disturbance  may 
be  determined  by  finding  the  stability  of  all  the  normal  modes  of  the 
chamber.  By  definition,  a  normal  mode  of  a  linear  system  varies  har¬ 
monically  in  time,  and  the  study  of  general  linear  stability  comes  down 
to  studying  the  behavior  of  harmonic  oscillations. 

Thus,  the  time  dependence  exp(iu>t)  is  assumed,  with  w  real  if  the 
oscillation  neither  grows  nor  decays  in  time.  In  order  to  display  more 
clearly  the  contributions  to  the  stability  of  motions,  it  is  more  con¬ 
venient  to  calculate  the  growth  or  decay  constant  t»,  defined  by  assuming 
the  complex  exponential  dependence  exp(  i  (ui-ia)t )  =  exp(iu>t)  exp(dft).  If 
a  is  positive,  the  motion  is  unstable;  processes  which  cause  a  to  be 
positive  are  said  to  be  destabilizing.  The  combination  k  =  (iu-ia)/a, 
where  a  is  the  speed  of  propagation  of  small  disturbances,  is  the  com¬ 
plex  wave  number.  All  analyses  of  the  linear  stability  for  combustion 
chambers  lead  eventually  to  results  for  k. 

The  general  analysis  is  based  on  the  equations  of  conservation 
lor  the  flow  in  the  chamber.  Although  the  procedure  is  a  bit  tedious, 
it  can  be  summarized  as  a  sequence  of  several  steps  which  are  impor¬ 
tant  to  recognize.  This  is  a  simple  way  of  emphasizing  the  important 
approx i mat  ions . 
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First,  in  a  more  general  vein,  the  formulation  of  the  problem  be¬ 
gins  with  the  notion  that  the  behavior  of  waves  in  a  combustion  chamber 
is  a  consequence  of  relatively  small  influence::.  The  idea  is  that  the 
structure  (spatial  distribution  and  time  dependence)  of  a  normal  mode 
for  the  real  problem,  including  combustion  and  flow,  is  "nearly"  that 
for  a  classical  acoustic  mode  in  a  chamber  having  the  same  geometry,  but 
with  no  combustion  and  flow.  If  the  frequency  of  the  nth  classical  mode 
is  u)  ,  then  the  constant  a,  and  the  change  of  frequency  u)-u>  for  the 
corresponding  mode  in  the  combustion  chamber,  are  small  compared  with 


or,  u»  -  u>  <<  u)  (21) 


It  is  not  assured  a  priori  that  this  inequality  is  true,  but  it  has  been 
verified  both  analytically  and  experimentally  for  most  situations  en¬ 
countered  in  practice. 

It  follows  from  this  approximation  and  the  assumption  of  linear 
behavior  that  there  are  two  small  quantities  characterizing  the  general 
problem:  the  amplitude  of  the  motions,  and  a  measure  of  the  perturba¬ 
tions  due  to  combustion  and  average  flow.  These  are,  respectively,  the 
Mach  number  M‘  for  the  amplitude  of  the  fluctuation,  and  a  Mach  number 
M  representative  of  the  average  flow  field.  Construction  of  a  suitable 
set  of  differential  equations  rests  on_a  formal  limit  process  defined  by 
the  relative  sizes  of  M'  and  M  .  If  M  =  0  and  M'  +  0,  all  influences 
of  combustion  vanish,  and  classical  linear  acoustics  is  recovered.  If 
M^  =  0,  and  M*o  remains  finite,  equations  for  nonlinear  classical  acous¬ 
tics  are  found. 

The  equations  for  linear  acoustics  with  combustion  and  flow  are 
defined  by  the  limit  process 

M*  ►  0,  M  1;  M’/M  *  0  (22) 
o  o  o  o 


This  implies  that  the  Mach  numbers  of  both  the  acoustic  and  average  flow 
fields  must  be  small,  but  that  for  the  acoustic  field  must  be  smaller. 
If  M  *  ‘  1  but  is  not  allowed  to  go  to  zero,  equal  ions  for  nonlinear 

acoustics  with  combustion  are  obtained.  Th»*  simplest  nonlinear  problems 
are  defined  by 

M‘  <<  1,  M  1,  M*  /M  ..  0(1)  (23) 

o  o  o  o 


This  limit  process  provides  the  basis  for  the  approximate  nonlinear 
analysis  described  in  the  following  section. 
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It  is  important  that  the  formal  basis  just  described  be  recognized, 
although  it  is  not  necessary  to  be  aware  of  all  the  details.  The  formu¬ 
lation  of  the  problem  eventually  determines  what  input  data  are  required 
for  numerical  computations. 

With  these  preliminary  remarks,  we  now  describe  briefly  the  pro¬ 
cedure  leading  to  a  formula  for  the  complex  wavenumber. 

1.  The  general  conservation  equations  are  written  first  for  the 
mixture  of  gas  and  particles.  Then  the  equations  are  split 
into  two  sets,  one  set  for  the  gas  phase  and  one  for  the  par¬ 
ticulate  phase.  This  step  defines  the  representations  of 
interactions  between  the  two  phases:  exchanges  of  mass, 
momentum  (the  drag  force),  and  energy  (heat  transfer). 

2.  Earlier  work  on  combustion  instability  was  based  essentially 
on  the  two  sets  of  equations  formed  in  step  1.  However,  for 
propellants  heavily  loaded  with  metal,  this  turns  out  to  be  a 
poor  strategy.  The  speed  of  sound  in  the  mixture  may  be  20% 
(or  more)  less  than  the  value  for  that  of  the  gas  only.  If 
the  two  sets  of  equations  for  the  two  phases  are  maintained 
separate,  the  change  in  the  speed  of  sound  is  computed  as  part 
of  the  analysis.  But  the  change  is  too  large  a  perturbation 
to  be  handled  accurately.  Thus  it  is  much  more  satisfactory 
to  incorporate  the  influence  of  particle  loading  in  the 
"zeroth  order"  problem  by  suitable  combination  of  the  two  sets 
of  equations;  this  procedure  is  described  in  Reference  6. 
The  result  is  that  the  speed  of  sound  in  all  that  follows  is 
r. he  speed  of  propagation  ^for  small  disturbances  in  the  gas/ 
particle  mixture,  denoted  a. 


J.  The  equations  are  then  linearized  according  to  the  process 
defined  by  equation  (22).  This  produces  a  set  of  first  order 
linear  partial  differential  equations  with  source  terms  pro¬ 
portional  to  the  Mach  number  of  the  average  flow.  (There  are 
some  contributions  associated  with  gas/particle  interactions 
which  do  nut  depend  on  the  average  flow,  hut  are  neve*  the  less 
small  perturbations.)  These  equations  can  be  combined  to  pro¬ 
duce  an  inhomogeneous  wave  equation  for  the  pressure  field, 
and  an  associated  inhomogeneous  boundary  condition  for  the 
normal  component  of  the  pressure  gradient 


vv 
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The  assumption  of  harmonic  time  dependence  is  now  introduced 
to  give  the  problem  to  be  solved, 


V^p  +  k^p  =  h 


(26) 


n  •  Vp  =  -f 


(27) 


At  this  stage,  the  wavenumber  k  is  unrestricted — it  is  not 
identified  with  any  particular  mode. 

A.  Equation  (26)  with  the  boundary  condition  (27)  can  be  solved 
by  use  of  a  Green's  function,  followed  by  an  iteration  pro¬ 
cedure.  However,  to  the  order  of  small  perturbations  (linear 
in  M  )  usually  treated,  the  following  simple  procedure  gives 
the  same  result.  The  perturbations  represented  by  fi  and  f  are 
necessarily  small  for  this  analysis  to  work  well.  Then  the 
mode  shape  p  and  wavenumber  k  cannot  be  very  different  from 
the  corresponding  quantities  (p  ,  k  )  for  the  classical  acous¬ 
tic  mode  in  a  closed  chamber  having  the  same  shape.  The  un¬ 
perturbed  problem  is  described  by  the  equations 


X2?  ♦  k2p  =0  (28) 

rn  irn  v 

n  •  Vpn  =  0  (29) 

In  a  certain  sense,  the  "difference”  between  the  perturbed  and 
unperturbed  problems  contains  the  information  we  seek  for  the 
gains  and  losses  ot  energy  associated  with  the  perturbations. 
This  is  contained  in  the  formula  for  the  wavenumber.  The 
"difference"  between  the  problems  is  formed  by  multiplying 

equation  (26)  by  p  ,  equation  (2R)  by  p,  subtracting  and  inte¬ 
grating  over  the  volume  to  give 

JIP^P  -  P^pJdV  ♦  (k2  -  k“)  /  P  p^tiv  =  jTipndV 

The  first  volume  integral  can  be  integrated  by  parts  by  use  of 
Green’s  theorem,  with  the  boundarv  conditions  (equation  (27) 
and  (29))  inserted,  leading  to: 

Jp^P  -  P^pJdV  =  f£)  |  r>(Vp  -  pvS,nl  •  n  dS  =  j>ndS 
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in  tne  second  volume  integral,  p  =  p  for  first  order  in  small 
quantities,  and  the  final  result  cannbe  written 


1 2,2. 
k  =  k  + 
n 


1 


»>2  ... 
p  dV 
rn 


dV+# 


f  Pnds} 


(30) 


This  result  is  the  basis  for  nearly  all  current  analyses  of 
linear  stability  in  solid  propellant  rockets  and  will  be  used 
here  to  examine  in  greater  detail  the  influence  of  linear 
velocity  coupling. 

A  A 

The  functions  h  and  f  are  defined  during  the  construction  of  the 
wave,  equation  (24),  and  its  boundary  condition.  They  contain  repre¬ 
sentations  of  all  processes  contributing  to  the  losses  and  gains  of 
acoustic  energy.  Only  one  part  will  be  examined  here.  It  arises  in  the 
boundary  condition  expressing  the  fact  that,  owing  to  the  response  of 
the  combustion  processes,  there  are  fluctuations  of  the  flow  inward  at 
the  surface.  The  origin  of  this  term  may  be  easily  demonstrated.  The 
linearized  momentum  equation  has  the  form: 


P 


0u; 

Bt 


V  p*  + 


For  harmonic  motions,  the  gradient  of  the  pressure  normal  to  the  surface 
is : 


n  *  ^p  =  -ipaku  •  n  (31) 


Actually,  there  is  another  contribution,  a  combination  of  part  h 
and  part  of  f,  which  must  also  be  included,  to  give  the  formula  for  k  : 


) 


-r> 


ndS 


(32) 


The  sign  has  been  changed  so  that  while  n  is  still  the  outward  normal, 
both  u  and  u  •  n  are  positive  inward.  it  is  often  useful  to  express 
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this  result  in  terms  of  the  mass  flow;  from  equation  (31)  of  Refer¬ 
ence  6, 


u  •  n  +  u  •  n  =  (m,  +  m,  )  (33) 
pi2  p  b  T 

The  particulate  loading  (p  / p),  the  mass  of  particles  per  volume  of 
chamber)  is  C  ,  and  AT  is  ^he  difference  between  the  fluctuations  of 
temperature  at  the  edge  of  the  combustion  zone  and  in  the  acoustic  field 
in  the  flow  outside  the  combustion  zone.  For  subsequent  discussion  we 
shall  ignore  the  temperature  differences  and  the  influence  of  particles, 
so  AT  =  =  0.  Hence,  equations  (32)  and  (33)  lead  to: 


(34) 


In  the  first  section,  equations  (12)  and  (13),  the  response  func-t 
tions  were  introduced  and  related  to  fluctuations  of  the  linear  burning 
rate.  What  appears  in  equation  (34)  is  the  fluctuation  of  mass  flux 
departing  the  combustion  zone  and  entering  the  chamber.  If  ont  accepts 
the  assumption  of  quasi-steady  behavior  in  the  gas  phase 


m/  =  p  r '  and 
b  *c 


%  =  pc? 


(35) 


where  p 


is  the  density  of  the  unburned  solid  propellant. 


Hence 


=  £ 

m,  r 
b 

and  the  two  response  functions  may  be  defined  as 


m.  /m,  m,  /  m,  m, /m, 

_  b  b  _  b _ b  .  p  _  b  b 

p'/p  p/p  V  (“'/a)v 


(36) 


Here  (u'/a)  represents  an  effective  velocity  fluctuation,  parallel  to 
the  surface,  which  approximately  accounts  for  rectification  and  the 
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threshold  velocity;  see  Figures  9  through  11  and  the  accompanying  dis¬ 
cussion  in  the  previous  section.  Because  the  analysis  described  here 
is  linear,  the  correct  nonlinear  influences  of  rectification  and  thresh¬ 
old  velocity  cannot  be  rigorously  accommodated  without  further  discus¬ 
sion.  These  nonlinear  effects  will  be  properly  treated  in  the  following 
section. 

For  the  sake  of  the  argument  here,  we  assume  that  pressure  and 
velocity  coupling  are  linearly  additive,  so  that  substitution  of  the 
definitions  of  equation  (37)  into  equation  (34)  gives: 


pa  JpJdV 


where  c^,  cv  are  real  constants  whose  values  are  unspecified.  They 
simply  express  the  relative  weighting  of  pressure  and  velocity  coupling. 
It  is  to  be  noted  that  p/p  in  equation  (37),  the  actual  pressure  fluc¬ 
tuation  has  been  replaced  by  (p  /p),  the  unperturbed  mode  shape.  This 
is  justified  in  the  previous  discussion  of  the  linearization  of  equa¬ 
tions  of  motion.  Because  m^/m^  is  first  order  in  small  quantities, 

retention  of  p/p  instead  of  PQ/p  in  equations  (37)  and  (38)  add  terms 

of  second  and  higher  order  which  must  be  dropped. 


Special  attention  must  be  paid  to  the  integrand  of  the  second  inte¬ 
gral  in  equation  (33).  The  circumflex  (A)  over  the  term  Rv(u'/a)v  means 

that_the  exponential  time  dependance  exp(iu)t)  must  be  removed.  Because 
(u'/a)v  is  nonlinear  when  rectification  is  accounted  for,  this  requires 

careful  interpretation  explained  below  following  equation  (54). 

Now  k  is  complex,  k  =  (u>  -  ia)/a,  and  its  square  is: 


k 


2 


a 


(39) 


The  bracketed  terms  in  equation  (38)  are  corrections  of  first  order  in 
the  Mach  number  of  the  average  flow.  Consequently,  both  a  and  the 
change  of  frequency,  u>  -  u)  ,  are  first  order,  and  again  because  higher 
order  terms  must  be  dropped,  equation  (39)  must  be  approximated  by: 


(40) 
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With  kR  =  wn/a>  the  real  and  imaginary  parts  of  equation  (38)  are: 


The  meanings  of  the  real  and  imaginary  parts  of  the  integrands  in  the 
second  integrals  are  explained  below. 


2.2  RAYLEIGH’S  CRITERION  DERIVED  FROM  LINEAR  STABILITY  ANALYSIS 

Before  proceeding  further  with  equation  (42),  it.  is  instructive 
to  examine  the  interpretation  of  a,  in  particular  its  relation  to 
"Rayleigh's  criterion."  According  to  the  definition  of  k,  the  depen¬ 
dence  of  the  pressure  on  time  is: 


p’  iakt  i(ui-ia)t  ut  iuit 

!r  =  e  =  e  -  e  e  (43) 

P 

Hence,  if  or  is  positive,  the  fluctuation  is  unstable  and  grows  exponen¬ 
tially  with  time.  The  time  averaged  acoustic  energy  is  propor- 

2 

tional  to  the  time  average  of  the  real  part  of  (p')  ,  so: 


2ot 

e 


Differentiate  with  respect  to  time  to  produce  the  useful  n*l.ition  for  u 


i  i  ■  £' 

“  ■  2  .It  <«) 

In  equation  (42),  tin*  denominator,  Jp“  dV,  is  proportional  to  -  £  ‘ 
and  the  numerator  therefore  represents  the  rale  of  increase  of  lime- 
averaged  acoustic  energy  in  the  chamber.  This  is  the  origin  of  the 
statement  that  tf  represents  the  difference  between  "gains*'  and  "losses" 
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of  acoustic  energy.  All  processes  contributing  to  the  stability,  or 
instability,  of  waves,  may  be  interpreted  in  terms  of  their  influences 
on  the  rate  of  change  of  acoustic  energy  in  the  chamber. 

The  criterion  formulated  by  Rayleigh  may  be  found  on  Page  226  of 
Reference  8. 

”If  heat  be  given  to  the  air  at  the  moment  of  greatest  conden¬ 
sation,  or  taken  from  it  at  the  moment  of  greatest  rarefac¬ 
tion,  the  vibration  is  encouraged.” 

The  idea  is  that  for  greatest  driving  of  the  waves — maximum  tendency 
towards  instability--heat  addition  must  occur  in  phase  with  the  pressure 
wave  and  in  the  vicinity  of  a  pressure  anti-node.  This  conclusion  may 
be  deduced  directly  from  equation  (30).  If  distributed  heat  addition, 
Q  (Energy/vol.-sec) ,  is  included  in  the  equation  for  conservation  of 
energy,  it  eventually  contributes  the  following  term  to  h’  (equa¬ 
tion  (46)  of  Reference  6): 


(h-)  = . 

^  ;heat  addition  -2'  y  /3t 

a  1 


(45) 


For  harmonic  motions, 


(fi) 


heat  add 


ition  *  “ 


(46) 


2 

According  to  equation  (30),  the  corresponding  part  of  k  is 


»2  =  k2  (>--■  ifc 


dV 


(47) 


The  real  and  imaginary  parts  are 


(48) 


(49) 
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For  a  given  mode,  a  is  increased  if  the  product  Q^Pn  is  increased. 

This  depends  only  on  the  real  part  of  Q,  the  part  of  the  fluctuating 

heat  addition  which  is  in  phase  with  the  pressure  fluctuation  and  Q  p 
is  obviously  greater  if  the  heat  addition  occurs  where  the  pressure  is 
greater,  i.e.,  at  the  anti-nodes. 


Thus,  Rayleigh’s  criterion  and  the  general  principle  cited  in  the 
first  section  are  both  contained  in  the  same  analysis.  These  are  merely 
two  ways  of  expressing  the  transfer  of  energy  required  to  drive  and  sus¬ 
tain  acoustic  waves. 


2.3  RECTIFICATION  AND  THE  EFFECTIVE  VELOCITY 

The  formula  (42)  is  the  general  representation  for  the  coupling 
between  surface  combustion  processes  and  acoustic  waves.  It  is  the  pre¬ 
cise  representation  of  the  discussion  given  in  the  introduction.  As 
argued  earlier  and  as  shown  in  equation  (42),  it  is  the  real  part  of  the 
response  function  for  pressure  coupling  which  determines  the  tendency 
for  instability  due  to  pressure  coupling.  The  second  term  in  equa¬ 
tion  (42)  contains  the  real  part  of  the  product  of  the  response  function 
for  velocity  coupling  and  the  effective  velocity. 

Because  of  the  way  in  which  the  analysis  is  constructed,  the  re¬ 
sponse  function  and  the  velocity  fluctuation  (u/a)  are  independent. 
Moreover,  according  to  equation  (9)  (u/a)  contains  only  a  part  which  is 
out  of  phase  with  p^/p.  Hence,  for  use  in  equation  (42), 


(50) 


as  shown  earlier  by  equation  (17)._  The  difference  between  (17)  and  (50) 
is  that  the  subscript  (  )  on  (u/a)  indicates  that  rectification  and  a 
threshold  velocity  are  to  be  accounted  for.  This  demonstrates  once 
again  that  for  velocity  coupling  it  is  the  imaginary  part  of  the  re¬ 
sponse  function  which  controls  the  tendency  for  instability. 

According  to  equation  (SO),  treatment  of  linear  stability  asso¬ 
ciated  with  velocity  coupling  comes  down  to  two  problems:  analysis  or 
measurement  of  the  imaginary  part  of  the  response  function;  and  determi¬ 
nation  of  the  effective  velocity  (u/a)v>  The  second  problem  is  essen¬ 
tially  a  matter  of  kinematics  which  has  been  treated  in  several  places, 
e.g.,  References  1  through  3  and  9  through  11.  Little  attention  will  be 
given  here  to  efforts  to  calculate  the  response  function  itself.  The 
intent  is  primarily  to  provide  the  basis  for  interpreting  data.  For 
this  purpose,  the  effective  velocity  must  be  known. 
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What  is  required  is  the  imaginary  part  of  (u/a)v,  that  part  which 
is  out  of  phase  with  the  pressure  fluctuation.  Consider  again  the  sim¬ 
plest  case  of  the  fundamental  longitudinal  mode  in  a  straight  chamber. 
Either  the  real  or  the  imaginary  part  of  the  complex  pressure  field  may 
be  used,  and  to  be  definite  we  use  the  real  part  of  equation  (2): 


=  e  cos  ^  cos  u)jt  (51) 

P 

The  corresponding  velocity,  parallel  to  the  lateral  surface  is  (5): 


u’  =  e  -  sin  ~  sin  u.t  (52) 

Y  l  i 


which  is  entirely  out  of  phase  with  the  pressure;  the  frequency  is  of 
course  =  ak  =  a/2L  for  the  fundamental  longitudinal  mode. 

The  effective  velocity  must  include  the  influence  of  rectification 
and  the  fact  that  these  may  be  a  threshold  velocity.  We  shall  list  here 
the  various  possible  cases,  but  the  details  will  be  worked  out  only  for 
the  case  when  the  threshold  velocity  vanishes. 

Three  velocities  are  important:  the  average  flow  speed,  u,  the 
acoustic  velocity,  u',  and  the  threshold  velocity  u  .  We  assume  that 
the  average  flow  is  from  left  to  right,  u  then  being  positive  always; 
the  threshold  velocity  is,  by  definition  positive.  In  steady  flow, 
erosive  burning  is  supposed  to  occur  whenever  the  average  speed,  u,  ex¬ 
ceeds  the  threshold  velocity.  For  unsteady  motions  we  assume  that  ero¬ 
sive  burning  occurs  whenever  the  resultant  velocity  u  ♦  u*  exceeds  u^ 
for  motion  to  the  right  and  when  -(u  ♦  u')  is  greater  than  u^  for  motion 
to  the  left.  This  becomes  clearer  from  the  diagrams  given  below. 

What  we  call  velocity  coupling  occurs  when  the  erosive  burning  at 
any  time  during  unsteady  motions  exceeds  the  erosive  burning  which  would 
take  place  under  the  average  conditions  existing  at  the  same  time.  The 
circumstances  when  there  is  no  threshold  velocity  are  most  readily 
understood.  There  are  then  two  possibilities:  the  magnitude  of  the 
acoustic  velocity,  |u'|,  is  either  greater  or  less  than  the  average 
velocity.  These  cases  are  sketched  in  Figure  12.  The  shaded  areas 
denote  the  velocity  effective  in  causing  erosive  burning  under  unsteady 
conditions.  Note  that  rectification  occurs  whenever  the  resultant  velo¬ 
city,  u  ♦  u',  becomes  negative,  or  the  total  effective  velocity  is 
|u  ♦  u'|. 
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mm 


V  / 


FIGURE  2.  The  Effective  Velocity  Due 
to  Rectification. 


Steady  erosive  burning  is  proportional  to  u  itself,  so  the  effec¬ 
tive  velocity  for  velocity  coupling,  when  there  is  no  threshold  velo¬ 
city,  is 


|u  ♦  u  '|  -  u) 


Before  considering  the  more  complicated  case  of  u  ^  0,  let  us  work  out 
the  consequences  of  our  discussion  so  far,  to  see  how  the  results  fit 
into  the  analysis  of  stability.  According  to  equation  (50)  and  earlier 
remarks,  we  need  that  part  of  equation  (53)  which  oscillates  out  of 
phase  with  the  pressure  fluctuation  hut  at  the  same  frequency.  To  find 
that  contribution,  the  formula  (52)  is  used  for  u'  and  (53)  is  expanded 
in  Fourier  series, 


(.  )  =  U'o  ♦  2J  Cn  cos  ( nuij  t )  ♦  21  sn  sintnujjt) 
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For  use  in  the  formulas  (41)  and  (42)  we  need  only  the  terms  which  os¬ 
cillate  with  frequency  u^: 


|  =  Cj  cos  u>jt  +  Sj  sin  uj^t 

v 


To  see  what  is  required  in  (41)  and  (42)  and  to  make  the  connection  with 
the  use  of  the  complex  notation,  write  this  as  the  real  part  of  the  com- 
lex  function 


iw,t  i(ui, t  -  n/ 2) 


(f)rc--  '  -v  1 


=  (C1  -  iSx)e 


iu>i  t 


Then  in  (41)  and  (42),  we  can  write 


Rv(!t)  =  (Rir)  +  i*in>cc,  -  iS,)e'Wlt 


=  (R 


(r) 


RU)S.) 
v  1 


RJr)S,)e 


iu.  t 


Hence  we  find 


(S5)a,  b 
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It  is  easy  to  show  that  the  result  (55)b  reduces  to  the  special 
case  (50)  when  u'  is  the  classical  acoustic  velocity.  The  complex  fora 
of  u'/a  is  then  given  by  (8), 


u*  .  e  .  •  5> 

—  -  Z  sin(nr)e 

a  Y  L 


and 


iiu.t 


3 

-  |  -  iR^e^ 


Hence,  corresponding  to  (55)a,  b,  we  have 


(f)l 


(r) 


=  R^  C-  sin (*?) 


The  first  of  these  is  exactly  (50)  because  -u/a^  =  (c/y)  sin(nz/L) 
from  equation  (9). 


Although  only  the  values  of  Aj  and  Bj  are  required  here,  we  shall 
evaluate  all  the  coefficients  in  (54).  The  foruaulas  for  the  coeffi¬ 
cients  are: 


C  = 
o 


d(u>jt) 


In 

Cn  =  In  f  (  d(u>j t )  (56)a ,b,c 

o  a  v 

n 

S  =  ~  fir  |  sin(nui.t)  d(u»,t) 
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The  difficult  part  of  evaluating  the  formulas  arises  from  the  term 
|u  ♦  u’|  in  (53),  so  this  will  first  be  treated  separately.  With 
(52),  this  factor  may  be  written 


7  (u  +  u* )  =  1  +  p  sin  u>.  t  (57) 

u 


where 


P  =  ^  sin  ¥ 
Yu 


(58) 


All  possible  cases  will  be  covered  by  supposes  first  that  p  >  1.  The 
function  (58)  appears  as  sketched  in  Figure  12(a)  and  expanded  in  Fig¬ 
ure  13  for  one  cycle,  0  <  u>t  <  In. 


KIGUKK  13.  One  Cycle  of  a  Rectified  Waveform  for 
Mean  Flow  in  the  Positive  Direction. 


It  is  clear  from  l he  sketch  that 

1 1  ♦  p  s i n  u>l  |  -  (  l  ♦  |i  si n  uj j l ) 


0  '  ujjt 


J  n 


- 1  1  ♦  p  sin  o^l) 


in  ,  in  A 

,t  -  €J  '  Wjt  *  ,t  ♦  <1 


(l  ♦  |i  sin  Ujl  ) 


In 


-  a  •  Uji  -  in 
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The  value  of  o  is  set  by  the  condition  1  ♦  p  sin  (ip  +  <0  =0,  giving 

cos  o  =  i  (p  >  1)  (59) 

For  p  <  1,  there  is  no  rectification,  and  o  =  0.  With  this  represen¬ 
tation,  the  following  integrals  may  be  found: 


2  n 

\nj  I1  +  P  ,in  “jtjdCuijt)  = 


-  2  c  +  fi 


I-50  +  5  tin  a 
71  71 


(60) 


271 


•j  j  |l  ♦  p  sin  uij 1 1 cos(nu)jt)d(u)jt) 


2  n  2 

-  -  cos  n  -  -  nn 

n  z  n 


~  sin  o|l  -  |  sin  oj  (n  =  1) 


(61) 


2n 


~  /l  +  P  sin  u)jt  sin(nu>jt)d(u»jt) 


i  1) 


4  2  1 

■z  sin  O  ♦  p  1  *  •  .  sin  2o|  (n  =  l) 

n  7i  n 


(62) 


Values  of  these  integrals  for  p  <  1  (no  reel i f ical ion)  follow  directly 
upon  setting  0=0. 
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The  integrals  over  u/a  in  (56)  contribute  only  to  C  ,  and  the  final 
results  for  no  threshold  velocity  (ufc  =  0)  are 

Cq  =  7 [l  -  ~  a  ♦  jjfi  sin  ol  -  u  =  ^  ~[p  sin  o  -  o)  (63) 

a  a 

Cj  s  0  (64) 


S 


ri 


sin  no  ♦ 


26  J  sin(n-l)a  sin(n+l)) 
n  \  (n-l)  (n+li  / 


(67) 


We  have  assumed  here  that  only  a  simple  sinusoidal  oscillation,  having 
frequency  u)j  is  present  in  the  chamber.  This  may  be  due  to  any  mode. 

Because  =  0,  the  driving  of  this  mode  itself,  due  to  velocity  coup¬ 
ling,  is  proportional  to  the  coefficient  S  ,  as  the  remarks  following 
(54)  having  shown.  For  use  in  (42),  we  therefore  write 


R 


v 


(r) 


sin  o  ♦ 


0  -  sin  2o 

JT 


>1 


(68) 


The  other  terms  S  cos(nujji)  in  the  series  expansion  (54)  represent  the 
influence  of  velocity  coupling  on  harmonics  of  the  mode  having  frequency 
u*.  Ve  shall  see  in  the  foi  lowing  section  how  these  are  incorporated  in 
an  approximate  nonlinear  analysis.  Appendix  A  contains  an  application 
to  T-burner  data. 

Now  suppose  one  is  concerned  with  the  general  ion  of  a  mode  which 
has  frequency  u  not  an  integral  multiple  of  the  !  requency  u*j .  For  ex¬ 
ample,  one  might  suppose  that  a  longitudinal  mode  may  excite  the  trans¬ 
verse  modes  _  in  a  cylindrical  chamber.  For  use  in  (42),  one  needs  that 

part  of  (u/a)  which  is  proportional  to  cos  w  t  with  u  t  u.  .  Because 
r  v  r  r  o  o  1 
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we  have  assused  steady  waves  throughout,  (54)  say  be  used  as  a  valid 
representation  for  all  tine,  -•  <  t  <  ♦«.  To  find  the  contribution  for 
any  frequency,  we  may  take  the  Fourier  transform  of  (54)  and  extract  the 
component  having  frequency  u>  .  Define  the  transform  pair 


(69) 


U(f)ei2nftdf 


(70) 


But  after  (57)  is  substituted  into  (65)  and  the  integrals  carried  out, 
we  find  that  U(f)  consists  of  an  infinite  series  of  delta  functions 
placed  at  f  =  iu^/271,  ±2  w  /2ft ,  ±3  id./2k,  . .  .*mu./2n.  There  is  no  com¬ 
ponent  at  the  frequency  u>  .  The  reason  for  this,  of  course,  is  that 
(u'/a)v  defined  by  (53)  an<f  expressed  as  (54)  is  periodic  with  frequency 
u>.  and  therefore  has  no  component  with  frequency  different  from  integral 
multiples  of  uij . 

The  important  conclusion  is:  an  unstable  longitudinal  mode  will 
excite,  due  to  velocity  coupling,  only  lontitudinal  modes.  The  non¬ 
linear  behavior  due  to  rectification  can  generate  only  integral  har¬ 
monies  of  the  primary  mode.  We  must  emphasise  that  the  only  nonlinearity 
accounted  for  here  is  rectification.  Nonlinear  behavior  associated  with 
the  flow  within  the  chamber  may  generate  other  modes,  but  the  matter  has 
not  been  studied. 

When  the  threshold  velocity  is  non-zero,  the  situation  is  consider¬ 
ably  more  complicated;  a  single  calculation  will  not  cover  all  possi¬ 
bilities.  Rather  than  include  the  calculations  here,  we  shall  merely 
sketch  the  cases  for  information  and  to  complete  this  discussion  of  the 
influences  of  rectification  and  threshold  velocity. 

It  is  convenient  to  classify  the  possibilities  into  two  groups, 
those  for  which  the  average  flow  speed  is  greater  than  threshold  aid 
those  for  which  the  speed  is  less  than  threshold.  First,  for  u  >  u( , 
there  are  three  cases  distinguished  by  the  magnitude  of  the  acoustic 
velocity.  These  are  sketched  in  Figure  14.  Again  the  velocity  effec¬ 
tive  in  causing  erosive  burning  under  unsteady  conditions  is  shown  by 
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the  shaded  region^.  Note  that  there  is  not  erosive  burning  if  the  re¬ 
sultant  velocity  u  +  u '  does  not  exceed  the  value  of  the  ^threshold  velo¬ 
city.  Steady  erosive  burning  is  now  proportional  to  u  -  ufc,  so  the 
effective  velocities  for  velocity  coupling  are  given  by  the  following 
formulas  corresponding  to  Figure  14(a),  (b)  and  (c): 

(a)  |u'|  <  u  -  ut 


[|u  +  u'l  -  utI  -  (u  -  ut) 


=  |u  +  u '  |  -  u  (71) 


(b)  (u  -  ut)  <  | u '  |  <  (u  +  ut) 


=0  (-ut  <  u  +  u'  <  ut)  (72) 

\  a 

V 

(f )  =  1 1  u  +  u'l  -  ut]  -  (G  -  ut) 

3 

V 

=  |u  +  u'l  -  u  (u  +  u'  >  u^)  (73) 
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FIGURE  14.  Rectified  Waveforms  When  the  Mean  Velocity 
Exceeds  the  Threshold  Velocity. 


NWC  TP  6363 


(c)  |  u " |  >  (u  +  ut) 


0 


(“Ut  <  u  +  u'  <  Ut) 


[|u  +  iT  |  -  ut]  -  (u  -  ut) 


(7'0 


=  I  u  +  u '  |  -  u 


(u  +  u'  >  ut)  (75) 


[|u  +  u'|  -  (u  +  ut)J  -  (u  -  ut) 


=  |  u  +  iT  |  -  2u 


(u  +  u'  <  -ut)  (76) 


The  second  group  of  possibilities  arises  when  the  mean-flow  speed  is 
less  than  the  threshold  velocity,  u  <  u  .  The  three  cases  are  sketched 
in  Figure  15.  Now  the  portions  of  each  Cycle  during  which  erosive  burn¬ 
ing  may  occur  are  much  smaller  than  for  the  previous  cases.  The  formu¬ 
las  for  the  effective  velocity  coupling  are 


(a)  |u'|  <  ut  -  u 


(M  =  0  (77) 

a  7 
v 

(b)  (ut  -  u)  <  u'  <  (u  +  ut) 


(ir)  =  o  (5  +  u'  <  ut)  (78) 

3 

V 


U  +  II  '  -  U 


(U  ♦  u'  >  U.) 


(79) 


FIGURE  15.  Rectified  Waveforms  When  the  Threshold 
Velocity  Exceeds  the  Mean  Velocity. 
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(c)  |u'|  >  (u  +  ufc) 


* 

M  =  o  (-ut  <  u  +  ii'  <  u  )  (80) 

a  / 

v 


(r~j  =  l«  +  u'|  -  ut  (u  +  u'  >  ut)  (81) 

3 

V 


(r~)  =  |  u  +  u'|  -  (u  +  ut)  (u  +  u'  <  -ut)  (82) 

'  a  ' 

v 


For  use  in  the  stability  analysis  the  effective  velocity  must  in 
these  cases  be  analyzed  into  its  components.  Only  the  part  proportional 
to  sin  u)t  is  required  to  determine  the  contribution  to  driving  the  fun¬ 
damental  mode.  As  an  unstable  wave  grows,  the  conditions  at  a  given 

position  in  the  chamber  (i.e.,  given  u  and  u^)  progresses  from  case  (a) 
in  Figures  14  and  15  to  case  (c).  And  of  course  at  a  given  location  in 
the  chamber,  the  relative  values  of  the  mean-flow  speed  and  threshold 
velocity  change.  Hence,  the  appropriate  formula  for  the  effective  velo¬ 
city  (equations  (71)  through  (82))  changes  with  time.  A  correct  analy¬ 
sis  of  stability  therefore  involves  a  considerable  amount  of  bookkeep¬ 
ing.  Some  relevant  numerical  results  have  been  given  in  Reference  5. 

It  is  apparent  from  the  preceding  remarks  that  because  of  rectifi¬ 
cation,  the  growth  constant,  a,  given  by  equation  (42)  with  (68)  used 
in  tne  second  term  becomes  a  function  of  the  amplitude  |u'|  or  |p'|  of 
the  oscillations.  Thus,  what  started  as,  and  is  strictly,  a  linear 
analysis  of  stability  has  apparently  become  a  nonlinear  analysis.  The 
qualification  "apparent"  should  be  emphasized,  however,  because  it  is 
obvious  from  all  of  our  previous  discussions  that  no  truly  nonlinear 
analysis  has  been  accomplished.  The  situation  is  better  viewed  at  this 
point  as  characterized  by  a  parameter  (here  the  response  function  for 
velocity  coupling)  which  changes  with  amplitude.  A  method  for  treating 
velocity  coupling  within  a  true  nonlinear  analysis  is  described  in  the 
next  section. 

In  any  case,  it  is  certainly  true  that  as  a  consequence  of  rectifi¬ 
cation  it  is  quite  possible  that  both  the  magnitude  and  the  sign  of  the 
growth  constant  change  with  amplitude.  A  change  of  sign  is  especially 
interesting  because  it  offers  the  opportunity  for  true  nonlinear  sta¬ 
bility.  It  is  possible  that  a  is  negative  at  low  amplitudes  but  posi¬ 
tive  at  high  amplitudes.  Under  these  circumstances  a  small  amplitude 
wave  is  stable,  but  a  large  amplitude  wave  is  unstable.  An  initially 
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small  disturbance  will  therefore  decay,  but  a  sufficiently  large  dis¬ 
turbance  will  grow.  This  is  an  example  of  truly  nonlinear  unstability, 
sometimes  called  triggering. 

Appeal  has  often  made  to  the  preceding  reasoning  to  explain  the 
dependence  of  instabilities  on  the  amplitudes  of  pulses  initiated  in  a 
motor.  While  the  argument,  or  at  least  the  essential  idea,  may  be  cor¬ 
rect,  there  is  no  proof  and  not  a  single  example  which  can  be  unambigu¬ 
ously  explained  on  this  basis.  Indeed,  there  is  in  fact  a  serious  lack 
of  calculations  which  can  be  viewed  simply  as  credible  estimates  of  the 
phenomenon.  There  is  an  opportunity  to  carry  out  some  useful  calcula¬ 
tions  to  provide  a  more  systematic  and  quantitative  basis  for  conclu¬ 
sions  drawn,  within  the  idealized  framework  described  here. 

To  treat  a  nonlinear  contribution  correctly,  it  should  be  placed 
within  a  nonlinear  analysis.  Some  "exact"  numerical  analysis  has  been 
done,  e.g.,  Reference  11  and  more  recently  Reference  12,  with  increasing 
success.  Observed  behavior  has  been  reproduced  qualitatively,  and  in 
some  cases  quantitatively,  for  nonlinear  phenomena  in  solid  rocket 
motors,  although  little  has  been  done  specifically  to  clarify  velocity 
coupling.  The  calculations  are  expensive  and  restricted  to  one  dimen¬ 
sional  problems.  Consequently,  approximate  analysis  continues  to  be 
very  attractive  for  many  applications.  One  approach  is  summarized  in 
the  next  section. 


3.  APPROXIMATE  NONLINEAR  ANALYSIS 

The  approach  described  here  has  been  reported  in  detail  in  Refer¬ 
ence  12,  although  no  attention  is  paid  there  to  velocity  coupling.  We 
shall  not  discuss  the  manner  in  which  the  approximations  are  carried 
out;  the  intent  is  only  to  note  that  a  useful  framework  is  available  for 
assessing  the  importance  of  nonlinear  velocity  coupling. 

Because  the  linear  analysis  described  above  has  been  widely  useful 
for  many  years,  the  nonlinear  analysis  has  been  constructed  to  be  a 
logical  and  clear  generalization.  The  first  step  is  to  approximate  the 
general  conservation  equations  by  treating  the  Mach  numbers  of  the 
average  and  unsteady  flows  as  small  parameters.  The  limit  process  was 
cited  earlier  in  equation  (23).  Power  series  expansions  of  the  depend¬ 
ent  variables  (pressure,  velocity,  etc.)  lead  to  a  system  of  equations 
which  can  be  manipulated  to  give  a  wave  equation  and  boundary  condition 
lor  the  pressure  fluctuation.  These  have  the  same  form  as  equations 
(24)  and  (25),  except  that  the  functions  h  anti  f  are  now  nonlinear  in 
both  the  pressure  anti  velocity  fluctuations: 


VV 


1  * 
_  i  ■) 

a"  a 


-  h 


ii  *  Vp  '  =  1 


(83) 

(84) 


I 
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The  pressure  and  velocity  fields  are  then  expanded  in  series  of  the 
normal  modes  of  the  chamber,  essentially  a  Fourier  analysis  in  space 
with  time-varying  coefficients.  After  substitution  of  the  expansion, 
equation  (83)  is  integrated  over  the  chamber  and  the  boundary  condi¬ 
tion  (54)  is  used.  The  result  is  a  set  of  coupled  oscillator  equations 
for  the  amplitudes  Hn(0  of  the  normal  modes: 

2 

+  “n  %  =  Fn  ^ 


where  is  the  coefficient  of  the  nth  term  in  the  representation 


The  forcing  function  is 


-2 

F  =  — — 

n  2 

pEn 


{/fipudV  +  (§fpndS} 


E  2  =  Jp  2dV  (88) 
n  Jrn 

When  h'  and  {'  are  linear,  the  equation  (85)  can  be  solved  by  assuming 
exponential  behavior  in  time,  and  the  formula  (30)  for  the  wave  number 
follows  directly. 

More  generally,  the  force  F  in  the  equation  for  the  nth  mode,  or 
oscillator,  depends  nonlinearly  on  rj  itself  and  also  the  amplitudes  of 
all  the  other  modes.  Thus  the  unsteady  motions  in  the  chamber  are  re¬ 
presented  as  a  system  of  oscillators  coupled  by  nonlinear  interactions. 
The  interactions  cause  a  flow  of  energy  between  the  oscillators  which 
represents  the  generation  of  harmonics  in  the  acoustic  field.  Note  that 
the  formulation  is  not  restricted  to  one-dimensional  problems  nor  to 
anv  particular  geometrical  conf igurat ion. 
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After  the  functions  are  set  by  specifying  the  functions  h  and  f, 

(85)  s  a  system  of  coupled  ordinary  differential  equations  which  can  be 
solve j  numerically.  It  is  not  particularly  cheap  to  do  so  because  the 
equations  are  second  order.  Considerable  simplification  is  achieved  by 
taking  advantage  of  the  fact  that  in  most  practical  situations  the  waves 
have  amplitudes  which  vary  quite  slowly  -  the  fractional  change  is  small 
in  one  cycle  of  the  oscillation.  A  convenient  procedure  to  follow  is 
based  on  the  method  of  averaging  which  reduces  the  system  of  second 
order  equations  to  a  set  of  first  order  equations,  two  for  each  member 
of  (84).  Write  the  amplitudes 


nn(t)  =  An(t)  cos(V  +  4>n(t)) 

=  A  (t)  sin  to  t  +  B  (t)  cos  u)  t 
n  n  n  n 


(89) 


where  A  ,  F  ,  A  ,  B  are  slowly  varying  functions  of  time  in  the  sense 
described  at5ove.n 

Eventually  the  first  order  equations  for  the  A  and  B  can  be  put 
In  the  form  n  n 


dA 
_ n 

dt 


a  A  + 
n  n 


0  B  +  g 
n  n  an 


(90(a)) 


dB 

n 

dt 


a  B 
n  n 


0  A  +  h 
n  n  n 


(90(b)) 


The  linear  portion  of  produces  the  terms  containing  and  0^  which 

are  now  strictly  constant  with  no  dependance  on  the  amplitude  of  oscil¬ 
lations.  All  nonlinear  behavior  is  contained  in  the  functions  g  and 

n 

h  .  There  is  no  difficulty  accommodat ing  nonlinear  velocity  coupling. 
No  calculations  based  on  this  method  have  been  published. 
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NOMENCLATURE 


a 

C 

m 

C 

n 

f 

h 

k 

k 

n 
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mb 

M 

o 

n 

h 

P 


6 

P 


6 

r 


H 

v 

ttS 


S 

n 

t 

T 


v 


Speed  of  sound 
Particle  loading 

Coefficient  in  Fourier  series,  equation  (54) 

Frequency  (cycles  per  second);  function  defined  by  equation  (25) 
Defined  by  equation  (24) 

Wave  number,  k  =  (u>  -  ia)/a 
Wave  number  of  the  n*"*1  normal  mode 
Length  of  a  chamber 

Mass  flux  of  material  leaving  a  burning  surface 
Reference  Mach  number 

Exponent  in  the  burning  rate  low,  r^pn 

Unit  normal  vector 

Pressure 

Small  change  of  pressure 

Linear  burning  rate 

Small  change  of  burning  rate 

Response  function  (or  pressure  coupling 

Response  function  for  velocity  coupling 

Element  ol  area 

Coefficient  in  Fourier  series,  equation  (54) 

Time 

Temperature 
Ye  loci ty 
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dV  Volume  element 

z  Coordinate  along  axis  of  a  chamber 

a  Imaginary  part  of  complex  wave  number 

p  Defined  by  equation  (57) 

Y  Ratio  of  specific  heats 

£  Amplitude  of  an  acoustic  pressure  oscillation 

K  Thermal  diffusivity  of  solid  propellant 

A  Wavelength 

p  Material  density 

u>  Frequency  (radians  per  second) 

-2 

0  Dimensionless  frequency,  Q  =  r  uj/k 

(  )  Average  value 

(  )'  Fluctuation  including  time  variation 

(  )  Fluctuation  without  time  variation 

(  Real  part 

(  ) ^  ^  Imaginary  part 

I  |  Absolute  value 
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Appendix  A 

REDUCTION  OF  T-BURNER  DATA  TO  OBTAIN  THE  RESPONSE 
FUNCTION  FOR  VELOCITY  COUPLING 

The  fluctuation  of  burning  rate  due  to  velocity  coupling  it  written 


(A-l) 


and  the  formula  for  the  growth  constant ,  the  second  part  of  equa- 
tion  (42)  with  c^  =  0  and  c^  =  1 ,  may  be  written 


a 


vc 


1 

2 


(A-2) 


The  circumflex  f~)  denotes  the  part  without  time  dependence,  and  super¬ 
script  r,  (  )K  ,  means  real  part,  that  part  which  is  in  phase  with  the 
pressure.  Thus,  the  pressure  fluctuation  having  amplitude  c  is  written, 
equation  (51): 


P  Kz 

=  C  cos  r-  COS  U)  t 

p  L  1 


(A-3) 


-  -  n  z 

p  =  cp  COS 


(A-4) 


and  p  means  p  at  the  frequency  iu  of  the  nth  normal  mode.  First  we 
examine  the  situation  when  rectif  ica*lion  is  ignored. 

If  u*  stands  only  for  the  fluctuating  velocity  associated  with  the 
acoustic  pressure  (A-3),  then  it  is  given  by  equation  (52), 


u ' 


„  a  ,nz. 
c  -  S.n(r) 


sin  u)^t 


(A-5) 


NWC  TP  6363 


which  is  90  degrees  out  of  phase  with  the  pressure  fluctuation.  The 
response  function  is  a  complex  quantity  which  may  be  written  according 
to  equation  (11)  as 


R 


v 


=  |Rv|[cos  t|)y  +  i  sin  <|)y] 


R 


V 


(A-6) 


In  complex  notation, 
p '  =  ep  cos 


(A-3)  and  (A-5)  are 

iu>- 1  iu>- 1 
kze  =  pe 


-  i(w  t-n/2)  iw.t 

u '  =  £  -  sin  kze  ~  ue 


(A-7) 

(A-8) 


where  u  is  pure  imaginary, 


u  =  -i£ 


a  . 

-  sm  kz 


(A-9) 


Thus  (A-3)  and  (A-4)  are  the  real  parts  of  (A-7)  and  (A-8).  The  wave 

number  k  is  complex,  k  =  (uj  -  ia)/a  of  the  imaginary  part,  the  growth 
constant  a,  is  given  by  the  formula  (A-2). 

For  use  in  (A-2),  with  u/a  entirely  out  of  phase  with  the  pressure, 
and  Rv  expressed  by  (A-6)  we  have  (50); 


(A-10) 


Comparison  of  (A-7)  and  (A-8)  gives  the  familiar  relation  between  acous¬ 
tic  pressure  and  velocity, 


u  - 


d£ 


PW. 


(A-ll) 


Substitution  of  (A-10)  and  (A- 12)  into  (A-2)  leads  to 


a 

vc 


-2- 

i a 

2  pwa 


dV 


dV 


(A-12) 


For  a  uniform  cylinder,  dV  =  Scdz  and  with  p 


p  ,  this  becomes 
rn’ 


(A-13) 


^2 

But  because  p'  has  the  same  value  at  both  ends  of  the  chamber,  the  inte¬ 
gral  vanishes: 


fd2i i 

dz 


dz  =  p2(L)  -  p2(o)  =  0 
*nv  *nv 


It  is  because  of  the  different  phases  between  the  velocity  and  pressure 
fluctuations  that  there  is  no  influence  of  velocity  coupling  in  this 
case;  see  Figures  5  and  7  for  examples. 

The  main  point  is  that  if  rectification  is  ignored  and  the  config¬ 
uration  is  symmetric  (this  was  implied  in  carrying  out  the  integral 
above)  then  there  is  no  influence  of  velocity  coupling.  So  now  we  must 
account  for  rectification. 

We  assume  that  the  threshold  velocity  is  zero  and  write  the  velo¬ 
city  effective  in  velocity  coupling  as  equation  (53): 

(  7" )  =  *  I  u  +  u '  |  -  ±  |  u  |  (A- 14) 

'  a  ■  a  a 

v 


The  first  term  |u  ♦  u'|,  accounts  for  the  rectification  with  unsteady 
flow,  and  the  second  term  ensures  that  the  effect  vanishes  when  there 
is  no  fluctuation  (u'  -  o).  We  must  now  determine  that  part  of  (A- 14) 
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which  is  out  of  phase  with  the  pressure  fluctuation,  that  is,  oscillates 
with  the  time  dependence  sin  tut,  and  has  the  same  frequency  as  the  pres¬ 
sure  fluctuation. 

To  accomplish  this  most  easily,  expand  (A-14)  in  Fourier  series, 
formally  written  as  equation  (54) , 


00  00 

(r~)  =  Cq  +  Cn  cos(nu)jt)  Sn  sin  niu^t 

a  v  n=l  n=l 

where  the  coefficients  are  given  by  the  formulas 


(A- 15) 


271  . 

co  =  kf(f)  (A‘16) 

v 


27t  - 

Cn  =  -  f  cos  (nujj t)  d(uj1 1)  (A- 17) 

J  a  v 

o 


.  £.1  » 

>  -  f  sin  (nuijt)  d(u)jt)  (A- 18) 

a  v 
o 

As  remarked  after  equation  (54),  only  the  coefficients  and  are 
required. 

The  coefficients  and  are  given  in  the  main  text,  equa¬ 
tions  (63)  through  (67)  for  the  case  when  the  threshold  velocity  is  as¬ 
sumed  to  vanish  and  the  mean-flow  speed  u  is  positive,  implying  flow  to 
the  right.  Equation  (66)  for  is 


S 


1 


u 


s{sin 


o*  + 


n 


sin  2 a 


ii 


a 


(A-19) 
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where  cos  o  =  1/p  when  p  >  1  and  a  =  0  if  p  <  1.  The  case  o  f  0  arises 
when  rectification  occurs,  otherwise,  p  <  1  according  to  the  definition 
(58), 


P  = 


£a 

Yu 


sin(j~) 


(A-20) 


That  is,  p,  is  the  ratio  of  the  local  amplitude_of  the  acoustic  velocity 
for  the  local  mean-flow  speed.  In  a  =  0,  =  up  and  (A- 10)  becomes 


R(i) 

_v_  d£  =  R(i) 


|  sin(Tl^) 


(A-21) 


Substitution  into  (A-2)  correctly  produces  (A- 13)  again. 

We  now  consider  the  procedure  for  calculating  a  associated  with 
oscillations  in  a  T-burner.  For  simplicity  we  treat  only  the  case  of  a 
burner  having  a  center  vent  and  symmetrical  configurations  of  grains. 
End  grains  may  be  installed  to  provide  average  flow  part  test  grains, 
and  additional  pressure-coupled  driving  of  instabilities.  The  following 
analysis  is  intended  only  to  give  results  for  the  contributions  due  to 
velocity  coupling  with  symmetrically  placed  cylindrical  test  grains.  It 
is  not  a  comprehensive  analysis  of  the  T-burner. 

For  each  test  sample,  the  distribution  of  mean-flow  speed  and  the 
acoustic  velocity  must  be  known.  The  latter  may  be  calculated  from  the 
formula  (A-4)  for  the  acoustic  pressure.  The  sketches  in  Figure  A-l 
show  typical  <iist ributions  of  the  acoustic  and  mean-flow  velocities  for 
samples  on  two  sides  of  a  T-burner.  The  acoustic  velocity  is  drawn  for 
the  portion  of  a  cycle  when  the  oscillatory  motion  is  everywhere  to  the 
right. 


Note  that  the  signs  must  be  correctly  accounted  for  so  we  must  con¬ 
sider  the  case  when  u  and  u*  have  opposite  signs.  For  u  negative,  re¬ 
plan*  ii  by  -  u  in  the  formula  (A- 14)  to  find 


1 


u 


u  II  =  I  p  sin  wt  -  1  I 


( A-22) 
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FIGURE  A-l .  Mean  and  Fluctuating  Velocities  for 
Two  Sample  Grains  in  a  T-Burner. 


where  p  is  defined  by  (60)  with  |u|  replacing  u.  A  sketch  of  (A-22)  for 
u  0  is  given  in  Figure  A-2.  This  corresponds  to  the  sketch  in  Fig¬ 
ure  13  for  u  positive,  and  is  obtained  by  reflecting  Ip  sin  io  t  +  1 1 
about  the  origin,  and  then  shifting  by  +2rt.  The  piecewise  representa- 


|p  sin  uit  -  l|  =  1  -  p  sin  u*1 1  0  <  w  t  <  a 


-  P  sin  iUjt  -  1  a  <  uijt  <  ^  +  o 


=  1  -  P  sin  aij t  ^  +  a  <  w,t  <  2n 

The  value  cf  o  is  the  same  as  for  the  case  u  >  0  because  here  the  con¬ 
dition  is  p  sin  (”  -  <j)  =  1  which  gives  cos  a  =  1/p  as  in  equation  (58). 
Again  the  interesting  cases  arise  for  p  >  1;  the  first  three  coeifi- 
c* fnLs  in  the  Fourier  expansion  are: 


u_<  0  j  |  >  1 


Cj=.  jj  IP  sin  o  -  o| 


NWC  TP  6363 


FIGURE  A-2.  One  Cycle  of  a  Rectified  Waveform  for  Mean 
Flow  in  the  Negative  Direction. 


C2  =  0  (A-24) 

Sj  =  -  |  ^  sin  o  +  p  |l  -  ^  a  -  ^  sin  2o|  |  (A-26) 

The  average  value,  Cq,  is  the  same  as  for  the  case  u  positive,  as  a  com¬ 
parison  of  Figures  13  and  A-2  shows  immediately.  However,  the  value  of 
Sj  tor  u  <  0,  which  is  the  case  in  the  right  half  of  a  T-burner,  is  ex¬ 
actly  the  negative  of  that  given  by  equation  (A- 19)  for  u  >  0.  Hence, 
for  use  in  (A-2),  the  formula  (55)  is 


a 


S  (left  half) 

S  (right  half) 


(A-2  7 ) 
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where 


S  =  -  sin  a  +  p<l  -  §  a  -  ^  sin  2a  > 
71  (  7t  71  ) 


(A-28) 


Strictly,  the  integral  in  (A-2)  must  be  carried  out  over  the  sample 
grains,  but  as  a  first  approximation  for  quick  reduction  of  data,  it 
is  adequate  to  use  average  values.  We_assume  that  the  samples  are 
at  z  =  L/4  and  z  =  3L/4.  By  symmetry,  u  has  the  same  values  at  the 
two  locations.  If  the  areas  of  the  samples  are  each  S  ,  then  the  sur¬ 
face  integral  in  (A-2)  is  approximately 


vr) 


-  l(r) 


p  dS  S  S.  RU'  (S.p  )  ♦  (S.p  ) 

lvv  ^  L/4  1  n  3L/4 


=  £pS.  u  S  cos  j  -  ju|S  cos 


^SlvV S 


With  S  given  by  (A-28) 


-  i(r)  p 

\(b)  =  t 

a  p 

v  J  r 


f  « *  4-r ' 


(A-29) 


2  2-2  ScL 

Because  Jp^  dV  =  £  p  y,  th*?  formula  (A-2)  is  finally 


S  u 


/?/  evw  aveW  b\  D(i)  4  .  A  J,  2 

V  =  '  Rv  RS'”0*  pf  -  n 

f  r  n  ' 


2 a  sin  2a 


(A- 30) 


Both  p,  equation  (A-20),  and  therefore  a  =  cos  l/p  are  based  on  the 
value  of  u  averaged  over  the  propellant  sample. 


S8 
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